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We study the r]' g*g* vertex for arbitrary gluon virtualities in the time-like and space-like regions, 
using the perturbative QCD hard scattering approach and an input wave-function of the ?7'-meson 
consistent with the measured rj'-y*-y transition form factor. The contribution of the gluonic content of 
the 77'-meson is taken into account, enhancing the form factor over the entire virtuality considered. 
However, data on the electromagnetic transition form factor of the r;'-meson is not sufficient to 
quantify the gluonic enhancement. We also study the effect of the transverse momenta of the 
partons in the rj'-meson on the r)' g*g* vertex, using the modified hard scattering approach based on 
Sudakov formalism, and contrast the results with the ones in the standard hard scattering approach 
in which such effects are neglected. Analytic expressions for the r] g* g* vertex are presented in 
limiting kinematic regions and parametrizations are given satisfying the QCD anomaly, for real 
gluons, and perturbative QCD behavior for large gluon virtualities, in both the time-like and space- 
like regions. Our results have implications for the inclusive decay B rj' X and exclusive decays, 
such as B —* ri'{K, K*), and in hadronic production processes + N{N) rj'X. 



PACS numbers: 14.40.Cs, 12.38.Bx, 12.39.Ki 



I. INTRODUCTION 



The vertex involving the coupling of two gluons and 
ry'-meson i^^/g.g. (g^, g|, m^,) (here, and g| represent 
the virtualities of the two gluons) enters in a number 
of production and decay processes. For example, the 
inclusive decays B Tj'Xs 0] and exclusive decays 
B ^ rj'K 1^, ^, involve, apart from the matrix elements 
of the four-quark operators, the transitions b ~f sg* , fol- 
lowed by g* rj'g |^, 6 ^ sgg followed by gg rj' 
as well as the transitions g*g* — > rj' and g*g — > rj 
Thus, a reliable determination of the vertex function 
(which can also be termed as the gluonic transition form 



v' ii, 
7' |. 



factor) Friigrg* {qi, g|, m^,) is an essential input in a quan- 
titative understanding of these and related decays. Apart 
from the mentioned i?-decays, the r]'g*g^*^ vertex plays 
a role in a large number of processes, among them the 
radiative decay J/tp — >■ ri'j and the hadronic production 
processes N+N{N) — > r/'+X, where iV is a nucleon. The 
QCD axial anomaly jsj, responsible for the bulk of the rj'- 
meson mass, normalizes the vertex function on the gluon 
mass-shell, yielding (0, 0, m'^,). The question that 

still remains concerns the determination of the vertex for 
arbitrary time- like and space-like virtualities, qf; i — 1,2. 
A related aspect is to understand the relation between 
the rj' g*g* vertex and the wave-function of the T^'-meson. 
Stated differently, issues such as the transverse momenta 
of the partons in the ?7'-meson and their impact on the 
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r]'g*g* vertex have to be studied quantitatively. 

While information on the rj' g*g* vertex is at present 
both indirect and scarce, its electromagnetic counterpart 
involving the coupling of two photons and the Ty'-meson, 
Frfi-yj, more generally the meson-photon transition form 
factor, has been the subject of intense theoretical and 
experimental activity. In particular, the hard scattering 
approach to transition form factors, developed by Brod- 
sky and Lepage , has been extensively used in studying 
perturbative QCD effects and in making detailed compar- 
ison with data A variation of the hard scattering 
approach, in which transverse degrees of freedom are in- 
cluded in the form of Sudakov effects in transition form 
factors [|ll|, |l2j, has also been employed in data analyses. 
It has been argued ^ that the Sudakov effects 
improve the applicability of perturbative QCD methods 
down to moderate values of and give a better account 
of data on meson-photon transitions [ p5[ . For a critical 
review and comparison of the standard (Brodsky-Lepage) 
and modified (mHSA) hard scattering approaches, see 
Refs. We note that either of these approaches 
combined with data constrains the input wave-function 
for the quark-antiquark part of the ?7'-meson. However, 
the gluonic part of the Ty'-meson wave-function is not di- 
rectly measured in these experiments and will be better 
constrained in future experiments sensitive to the rj' g*g* 
vertex. 

A closely related issue is that of the rj ~ rj' mixing. 
There exist good theoretical and phenomenologi- 
cal ijl^ reasons to suspect that the frequently assumed 
pattern of the decay constants /p, defined by the equa- 
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tions 



{0\jl,\P{p))=ifl,p^, (/c = 8,l; P = 11,11'), 

(1.1) 

where J^^ and J^^g denotes the SU(3)i^ octet and singlet 
axial- vector currents, respectively, do not follow the pat- 
tern of state mixing. Defining the decay constants as fll] 



= /S COS 6*8, 

= /s sin 6*8, 



-fi sin 6*1, 
/i cos 6*1, 



(1.2) 



the angles di and (?8 are found to differ considerably due 
to non-neghgible SU(3)F-breaking effects |l^, |9|. In 
contrast, it is natural to expect that the state mixing 
is determined essentially by a single angle, as the state 
ric is far too heavy to be significant in the state-mixing in 
the ri — i]' complex. A growing consensus is now emerg ing 
in favor of the two mixing-angle scheme of Lcutwyler |lq| 
for the decay constants in Eq. (1.2) |l^, To be pre- 
cise, we shall use the state-mixing scheme of Feldmann, 
Kroll and Stech |2^, which is consistent with the two 
mixing-angle scheme of Ref. in the current-mixing 
basis. 

Based on the foregoing discussion, it is natural to use 
the hard scattering approach to study the ri' g* g* vertex, 
incorporating the information on the wave-function and 
the mixing parameters entering in the ?/— ry' complex from 
existing data involving the electromagnetic transitions. 
A first step in this direction was undertaken recently by 
Muta and Yang They derived the ri'g*g vertex in 

the time-like region in terms of the quark-antiquark and 
gluonic parts of the 77'-meson wave-function, taking into 
account the evolution equations obeyed by these par- 
tonic components |^. In this paper, we also address 
the same issue along very similar lines. We first red- 
erive the ii'g*g vertex, pointing out the agreement and 
differences between our results and the ones in Ref. . 
The latter have to do with the derivation of the lead- 
ing order perturbative contribution to the gluonic part 
of the ii'g*g vertex, and the use by Muta and Yang ||2^ 
of the anomalous dimensions derived in Ref. [ p2| in the 
evolution of the wave-functions. For the anomalous di- 
mensions, we use now, correcting a similar mistake in 
the earlier version of this paper, the results derived in 
Refs. [E3l_E4|, which are at variance with the ones given 
in Ref. p2f7but which have been recently confirmed by 
Bclitsky and Miiller [^. Making use of the ri'j*j data 
to constrain the ry'-meson wave-function parameters, we 
find that the gluonic contribution in the yy'-meson is very 
significant. We then extend our analysis to the case when 
both the gluons are virtual, having either the time-like 
or space-like virtualities. We also present a number of 
results on the vertex function Fr^' g* {qf , q2 , m'^,) in the 
asymptotic region, which have not been presented ear- 
lier to the best of our knowledge, and are of use in fu- 
ture theoretical analysis, in particular in the B — > ri'X 
and B ri'{K,K*) transitions. We study the effects 
of the transverse momentum distribution involving the 



constituents of the ry'-meson and take into account soft- 
gluon emission from the constituent partons by including 
the QCD Sudakov factor, following techniques en vogue 
in studies of the electromagnetic and transition form fac- 
tors of the mesons ||l^, |l^, |2^. We show the improve- 
ments in the il'g*g^*^ vertex function due to the inclusion 
of the transverse- momentum and Sudakov effects, which 
are particularly marked in the space-like region, improv- 
ing the applicability of the hard scattering approach to 
low values of Q^. These effects have a bearing on the 
hard scattering approach to exclusive non-leptonic de- 
cays the importance of the transverse-momentum 
and Sudakov effects in the decays B ^ tttt and B Ktt 
has also been recently emphasized in Ref. Finally, 
we also derive approximate formulae for the ii'g*g^*^ ver- 
tex, which satisfy the axial- vector anomaly result for on- 
shell gluons and the asymptotic behavior in the large-Q^ 
domain, determined by perturbative QCD. 

This paper is organized as follows: In section ||, we 
specify the 77 — 77' mixing formalism and the evolution 
equations for the qq and gluonic wave-functions of the 
77'-meson. In section ID, the qq contribution to the off- 
shell ii'g*g* vertex is worked out in the hard scatter- 
ing approach. The corresponding contributions from the 
gluons are presented in section IV . In section ^ we im- 
plement the transverse-momentum effects and derive the 
Sudakov-improved il' Q* g* vertex. Numerical results are 
presented in section VI, where we al so show comparison 
with existing results. In section VII, we give simple for- 
mulae for the i]'g*g* vertex, which interpolate between 
the well-known lim iting cases. Our main results are sum- 
marized in section VIII. Appendix ^ contains the solu- 
tion of the evolution equations for the wave functions 
<j)'^'^\x,Q) and <j)'^^\x,Q), and the function J{uj,ri) intro- 
duced in the derivation of the ri' g* g* vertex in section III 
is given in Appendix 



II. 



-MESON WAVE-FUNCTION 



The 7y'-meson is not a flavor-octet meson state and, 
hence, in addition to the usual quark content the 77'- 
meson wave-function has a gluonic admixture. In prin- 
ciple, there also exist a cc component in the 77'-mcson 
wave- function, but it has been estimated to be rather 
small in well-founded theoretical frameworks |2^, p9[ , 
and hence ignored. We take the parton Fock-state de- 
composition of the 77'-meson wave-function as follows: 



I77' > = sin0|77^ > +cos(j)\ii'^ >+K >' 



(2.1) 



where the SU(3)i;' symmetry among the light u, d and 



s quarks is assumed so that [77' > > 



dd > /V2 



and \rig > ^ \ss > are the quark Fock states and (f> is 
the mixing angle; |r/^ > ~ \gg > is the two-gluon Fock 
state. In Refs. [^2|, ^ ^ the eigenfunctions of the mix- 
ing quark \qq > and gluonic \gg > state: 
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(2.2) 
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V2f, 



fs COS ( 



(2.3) 



where fq and fs are the decay constants of \ri'^ > 
and 1 77'^ >, have been calculated by solving the evo- 
lution equations. The result for the quark and gluonic 
components is presented as infinite series of the Gegen- 
bauer polynomials of the indices 3/2 and 5/2 |59 and 



can be found in Appendix Previous analyses show 
that it is a good approximation to consider only a few 
first terms in the expansion of the quark and gluonic 
wave- function. Here, we shall keep the leading two terms 
in the expansion for cf)^'^') {x , Q) , and keep only the first 
term for (l)'^^\x,Q): 



^(«)(a;,Q) = 6xx<^ 1 + 



(9) 



{x, Q) = xx{x — x) 



(9) f as[Q^) 



B 



(9) 



15 Vas(A*o) 
,(9) f "s(<3^) 



5B' 



(1 — 5xx) 



(2.4) 
(2.5) 



where x and x = 1 — x are the energy fractions of 
two partons inside the yy'-meson, > is the en- 
ergy scale parameter, and ~ 0.5 GeV is the typical 
hadronic energy scale below which no perturbative evo- 
lution takes place. It is seen that in the limit ^ 00 



the quark wave-function (2.4) turns to its asymptotic 
form (/)as (a;) = 6xx (the same asymptotic behavior as the 
pion wave- function |^ due to i ts quark content), while 
the gluonic wave- function ( |2.5| ) vanishes in this limit, 
</>as^(a;) = 0. The coeffici ents of the expansion of the 
wave- functions (2.4) and ( ^.5D are calculated by using 
perturbation theory and include the effective QCD cou- 
pling as{Q^), which in the next-to-leading logarithmic 
approximation is given by pl| 



47r 



/3oln(Q7A2 



1 



2/3i lnln(QVA^; 
/32 ln(gVA2) 



(2.6) 



where /3o = 11 - 2n//3, /3i = 51 - 19n//3, A = Aqcd 
is the QCD scale parameter, and f^/ is the number of 
quarks with masses less than the energy scale Q. In the 
energy region fhc < Q < ffibj where ific = 1-3 ± 0.3 GeV 
and ifib = 4.3 ± 0.2 GeV are the charm and bottom MS 
quark masses, respectively, we will use the value of the 



dimensional parameter A^l^ 

MS 



280 MeV corresponding 



to four active quark flavors The requirement that 

the effective QCD coupling is continuous at the flavor 
thresholds leads to a change of A in the threshold energy 
regions where nf changes. This implies that for Q < rhc, 

the value of A is rescaled, yielding Ai-L — 333 MeV. 



III. QUARK-ANTIQUARK CONTRIBUTION 

TO THE v'9*9* VERTEX 

The diagrams depicting the quark-antiquark content 
of the 'rfn*g* vertex (or transition amplitude) are shown 
in Fig. n. The invariant amplitude corresponding to the 
quark contribution to the rj' g*g* vertex in the momentum 
space can be defined as: 



M 



(9) _ 



C 
4iV. 



dx^^''Hx,Q)Sp[j5{p-m^,)T^^>l (3.1) 



(9)1 



where and m,,/ are the four-momentum and the mass 
of the ry'-meson, respectively, p = Ptil^i Iti and 75 are 
the Dirac 7-matrices, and is the hard amplitude con- 
nected with the effective quark-antiquark-two-gluon ver- 
tex by the following relation: 




P^l^ = ^o^n y^r\^P. ^P. -91,-92) el* el 



(3.2) 



where the summation over the quarks' colors a and /3 
allows to get a color-singlet meson state and x and x 
have been defined earlier in Sec. In the leading order 
in the strong coupling a^, the effective quark-antiquark- 
gluon-gluon vertex is presented in Fig. ^ and has the 
form: 



FIG. 1; Leading Feynman diagrams contributing to the quark 
part of the r] g* g* vertex. 
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ii - 42)9^il^ + {qi + 292)^71^ - (2?! + 92)1/ 7m] 



(3.3) 
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FIG. 2: The effective quark-antiquark-gluon-gluon vertex and 
leading order contributions. 



{tbta)f3a 



(, , ^ IpXh + ^2 + mq)j^ 
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where ta {a — 1, . . . , — 1) are the generators of the 
color SU(iVf,) group, fate a-re its antisymmetric structure 
constants, g^^ = diag(l, —1, —1, —1) is the metric tensor, 
rriq is the quark mass, and ^ = 9/^7^ • In the calculation of 
the amplitude, we have neglected all the quark masses niq 
as they are much smaller than the energy scale parame- 
ter Q. 

Defining the quark contribution to the rj' g*g* vertex 

f''^ as- 
^n'g'g* as. 

(3.4) 

the result of the calculation is: 



1 



xq\ + xq2 — xxm^, + ie 



- — h (a; x) 



(3.5) 



Note that the quark wave-function of the r/'-meson satis- 
fies the symmetry condition 0'^'^^ (x, Q) = <j)^'^\x,Q) |2|, 
[24II . When one of the gluons is on the mass shell, 
for example — 0> the quark part of the t]' g*g* ver- 
tex, F^.g, (9i, 0, m^,), agrees with the one presented in 
Ref. g 

In the case when both virtualities of the gluons are of 
the same signs (space-like < 0, < or time-like 
q\ > 0, (7I > 0) the quark contribution to the rj' g*g* 



vertex can be presented in the form: 
47ra,(Q2) 2C 



(3.6) 



(1 - 2xxT]) <j)^'i\x,Q)dx 
(1 — 2xxrj)'^ — tjj'^{x — xY + 



where q^ = if + Q2 ^^e total gluon virtuality, u = 
(ql — q^)lq^ is the asymmetry parameter having the val- 
ues in the domain — 1 < w < 1, and r\ = m^i/q^ is the 
scaled Ty'-meson mass squared. As a typical mass scale Q 
which defines the value of the strong coupling as{Q^), 
it is natural to consider the absolute value of the total 



The quark contribution, F^'^l*g* 



virtuality — \q^ 

corresponding to k eeping the first two terms in the quark 
wave- function (2^) is: 



Fy,>g. {Q\^, V) = F^'^ {Q',u:, r,) + F^"^ {Q\ ,7) 
47ra,(Q2) QC ' 



6B. 



(q) ( a,{Q^) 



(3.7) 



(9) 
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where the functions Gq'"* and G^'' are: 



1^ , 

-IH In 

2r] 



1+UJ 



l-OJ 



UJ 



21 



G2 i^,V) = -r+7, 2" 



UJ 

2^ 







In 


l+LO 






1~UJ 




77^ 





jj3 



(3.8) 
(3.9) 



J(w,77). 



Notice that G'f\uj,ri) {i — 0,2) are symmetric in 
their first argument under the interchange u — s- —uj: 



G'f\—uj,ri) = Gf'{uj,r]). The function J(aj, r/) is defined 
and analyzed in Appendix 

In the space-like region of the gluon virtualities, where 
77 < 0, the rj' g* g* vertex is real as in this case the quark 
propagator does not have poles which is clearly seen 
from Eq. (B.6). In the time- like region of the virtuali- 
ties, for those values of the asymmetry uj and the rel- 
ative Ty'-meson mass squared rj for which the inequality 
w^ + (l — 77)^ > 1 is satisfied (Fig. an imaginary part in 
the vertex function appears due to the function J{uj,r]) 
as the result of the ie prescription of the quark propa- 
gator (see Eq. (|3.5[)). Notice that if one of the gluons is 
on the mass shell, for example the second one, = 0, 
the imaginary part is nonzero for any value of the other 
virtuality in the region > 7ti^, . For large gluon vir- 
tualities the imaginary part of the vertex becomes large 
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FIG. 4: The functions fo{i-^), /2(<^), and (?2(<j-') describing the 
large asymptotics of the ri' g* g* vertex, with the two gluon 
virtualities having same signs. 



FIG. 3: The uo ~ -q parameter plane of the rj'g*g* vertex. 
The regions over which the vertex has an imaginary part are 
indicated by the gray color. 



due to the linear dependence on of the leading contri- 
bution, Im Fq'^'' ~ , and the cubic dependence of the 

next-to-leading contribution, Im .Fj'''' ^ {q'^)^- This rapid 
increase of the imaginary part of the rj' g* g* vertex func- 
tion with the total gluon virtuality seems unphysical and 
the natural prescription in the Brodsky-Lepage approach 
is to drop the imaginary part. The physical interpreta- 
tion of this procedure will be discussed in Sec. 0. In the 
following, we shall drop the imaginary part of the rj' g* g* 
vertex in the Brodsky-Lepage approach. 

It is useful to find the asymptotics of the quark contri- 
bution to the r]'g*g* vertex at large {\r]\ ^ 1 or equiv- 
alently ^ m^, ) for arbitrary values of the asymmetry 
parameter lo. As in the case of the 7r7*7* transition form 
factor the function describing the i]' g*g* vertex decreases 
like - l/g^: 



F, 



9) 



(Q^ 0) ~ F^''^ {Q\uj, 0) + F'^'" (Q^ uj, 0) 



6B. 



(3.10) 
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(9) 
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troduced above are defined as follows: 



/o(^) = — 



2uo 



In 



l-LU 



(3.11) 



The functions fo{uj) and /2(t^) arc displayed in Fig. ^. As 
both these functions are symmetric in their arguments: 
fi{—uj) — ,fi{uj) {i — 0,2) we present them in the re- 
gion < w < 1 only. These functions take their maxi- 
mum values at the extremal values of the asymmetry lu: 
/o(±l) = 1 and /2(±1) = 1/6. At small values of the 
asymmetry parameter (|w| ^ 1), /2(tf) has a quadratic 
behavior /2('jj) — 4a;^/105 and goes to zero with iv, while 
fo{uj) is finite at a; = 0: /o(0) = 2/3. This implies that if 
the gluon virtualities are comparable to each other, i.e. 
qf ~ q2, the next-to-leading correction has an additional 
suppression given by the function f2{^)- 

When the virtualities of the gluons have opposite signs, 
the typical scale parameter can be defined as Q'^ — \q\~ 
In this case the quark contributions to the rj'g*g* 
vertex is: 



47ra,(Q2) 6C 1 



Nr. UJ 



\UJ OJ 



(3.12) 
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(9) f asiQ^) 
as(/io) 



B. 
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where in this case now q^ 



92^ 



w = (9? + g|)/9' 



where the two functions of the asymmetry parameter in- is the asymmetry parameter, the relative Ty'-meson mass 



6 



squared is 77 = mf , /q'^ , and the functions G^"^^ are defined 



by Eqs. and (^. 

In the region of large the quark contrib ution to the 
q* vertex has the same 

1/Q2 behavior ( pO| ) as m 
the case of the gluon virtuahties of the same signs, but it 
has a different dependence on the asymmetry parameter 
fi{uj) fi{uj) (i = 0, 2) where: 



The function /o(i^) is defined in Eq. (3.11), and the curves 
corresponding to these functions are presented in Fig. |^. 
Note that these functions are also symmetric in the asym- 
metry parameter: fi(—uj) = fi{'-^), like the fmictions 
encountered earlier. 



foico) 

/2H 



= 2 - C^Vo(^), 



(3.13) 



= -^ + ^ + ^(l-5c.Vo(-). 



J 



If one of the gluons is on the mass shell, for example 
the second one (ql = 0), the leading and next-to- leading 
contributions are simplified: 



F(«)(g?,0,m^,) = 
F^'^\qlO,ml,) 



-^--P^|l + 4lnfl-4^ 
Nc mf, \ qf 



'V' 

47ra,{\qf\) 6C 

Nr. 



m 



n 



a 



X < - 



1 bql , b{qlf , ql 



as 
,2 



48 I \ 



15 



6 2to, 
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rat. 



m'. 
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101 

81 
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(3.14) 
(3.15) 



In the limit of the large gluon virtuality ( | | 3> mfj, 

^^'r r («i ' 0' 0) ^2 ]^ 



,,) the quark contribution to the rj' g* g* vertex has the form: 



1 + B^^^ 



fas{\ql\)\ 


48 / % 


(asm)\ 
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101 
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When the gluon 



= rat 



virtuality is time-like q\ > 0, the quark contribution has a logarithmic divergence near 



n'9'9* 



ql-rn.1 



-l!^^^^|l + lnfl-4 



(3.16) 

the threshold 

(3.17) 
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For the space-lik e glu on vir tuali ty < it is worth 



noting that Eqs. ( 3.14 ) and ( 3. IE ) are suitable when the 
gluon four- momentum satisfies the condition: l^i | > rn^y' 
In the case of the smaller absolute values of the gluon 
virtuality the jy'-meson mass, m^/, becomes the largest 
scale parameter. The strong coupling as should be esti- 
mated at that value of the scale parameter which corre- 
sponds to the following change in Eqs. (3.14) and (3.15): 



IV. GLUON CONTRIBUTION TO THE r]' g* g* 
VERTEX 



The invariant amplitude of the gluon contribution to 
the ri'g*g* vertex can be defined as: 



iC 



dx0(^')(x,Q)£"'3/^"^gi„g2/3 



(9) 



(4.1) 

where T-^^ is the hard amplitude of the gluonic content 
shown in Fig. In the Brodsky-Lepage approach it is 
connected with the effective four-gluon vertex by the fol- 
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FIG. 5: The functions fo{i^), /2(i^), and g2{i^) describing the 
asymptotics of the rj'g'g* vertex function at large Q^, with 
the two gluon virtualities having opposite signs. 





FIG. 6: Leading order contribution to the gluonic part of the 
v'g*g* vertex. 



lowing relation: 



pa 



(4.2) 



where the summation over the gluons' colors c and d 
allows to get the two-gluon color-singlet state. In the 
leading order in the strong coupling as the effective four- 
gluon vertex has the contributions from the diagrams of 
the four-gluon annihilation shown in Fig. ^ and can be 
written in the form: 



F„X5(9l,'?2,g3,94) 



47ras fahnfcdn 



{{qi 



q2) 



(4.3) 



{qi + 92)^ + ie 

X [9a-fgi3s - 9as9p-i\ + Qafigis ((^i - q2){qs. - qi)) 
+gaf3 [(93 + 294)7(91 - q2)s - (2g3 + 94)5(91 - 92)7] 

+575 [{qi + 292)0(93 - 94)/3 - (291 + 92)/3(93 - 94)0] 

-1-507(291 -t- 92)/3(293 + 94)<5 + g0s{qi + 292)0(93 + 294)7 

-505(291 + 92)/3(93 + 294)7 - 5/37(91 + 292)0(293 + 94)5} 

+ {(a, a, 91) <-> (c, 7, 93)} + {{a, a, 91) <-> {d, S, 94)} . 
Let us define the gluonic contribution, F^y ^.^, , to the 




FIG. 7: The effective four-gluon vertex and leading order con- 
tributions. 



v' 9* 9* vertex as: 



M^"' = -T-^iYg'g' Oab£ £p £y 9ip92(T- 



(4.4) 



The result of the calculation is: 



47ras(g2) C 



- / dT(/.(f)(x,Q) 



xql 



xql 



(1 -I- xx)m?^y 



xql 



xql 



xxmz 



{x 



(4.5) 



Note that the gluonic wave-function of the jy'-meson 
satisfies the antisymmetry condition ^'^^ (x, Q) = 
-(j)'^3)(^x,Q) m, H, Hi. It implies that if the relative 
sign in the brackets were a (as given in Eq. (6) in 
ReL H^) the gluonic contribution to the rj'g*g* vertex 
would vanish identically. We disagree with the result 
stated in ReL ||2^ on this point. 

When the gluon virtualities are of the same signs, the 
gluonic contribution to the f]' 9*9* vertex can be pre- 
sented in the form: 



{Q^^,v) = 



47ras(Q^) 



2Cw 



(4.6) 



(x - 5) [1 - (1 -h 2xx)^\ 0(f) (x, Q) dx 
(1 — 2xxrj)'^ — cj2(a; — x)^ + ie ' 



where the total gluon virtuality 9^, the asymmetry pa- 
rameter cj, and 77 are the same as the ones entering in 



Eq. (3_^) and we have set the scale — |9^|. After the 



integratio n, a nd taking into account the gluonic wave- 
function (2.5), the result is: 



f(?].^.(Q^c.,,7) = 
(q) fasiQ'^ 



'-^%Gi^\.,,) (4.7) 



16B 



as(Mo) 



5B}, 
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where the function Gj^^ has the following expression: 



r^(Q), N I 5 2 4cj2 

G\'\u,ri) = uj{- + ^ 

' 3 r/ 77^ 



(4.8) 



r 21 




4c^2- 


In 










l-LO 






if' 





8uf 



„4 



J(tj,7;) 



Note that this function is antisymmetric in its first ar- 
gument under the change to — 



As in the case of the quark-antiquark contribution to 
the rj' g*g* vertex, the gluonic contribution also contains 
an imaginary part due to the function J(ti;, 77) as a result 
of the ie prescription of the gluon propagator. The imag- 
inary part of the gluonic contribution has a cubic depen- 
dence on the total gluon virtuahty, ImF^fj.^. ~ {(fY, 
in complete analogy with the next-to-leading quark con- 
tribution. We use the Brodsky-Lepage prescription dis- 
cussed in the preceding section and drop the imaginary 
part in the analysis of the vertex function. 

The gluonic contribution to the t]' g* g* vertex has the 
same as ympt otics at large {rj 0) as the quark 
contribution (3.1C): 



(4.9) 



IQB. 



(g) / as{Q^) 



3/0 (c^) - 2 



+ 5M 



as(/i§) 



(4.10) 



where the function /o(w) is defined by Eq. ( ^.11 ). 
Note that the gluonic function (72 (w) is antisymmetric: 
52 (—w) — — g2(^), in contrast with the quark functions 
fi{—Lo) = fi{i-o) {i = 0,2). The curve corresponding to 
52 (w) is presented in Fig. ^. It takes its maximum and 
minimum values at the borders of the argument domain: 
52 (±1) = ±1/6 and 52(0) — 0, respectively, because of 
the antisymmetry condition. Again, as in the next-to- 
leading order quark contribution to the rj' g*g* vertex, 
the gluonic contribution has an additional suppression 
due to the function 52 (w) for the gluons with comparable 
virtualities. 

For the case when the gluons virtualities have opposite 
signs {q\ > 0, q| < or ql < 0, g| > 0), ihejl^g* vertex 



can be presented in the form given in Eq. (4.7) but with 
a different dependence on the parameters w and rj: 



Gi'\io,rj)^G 



(9) 



1 



(4.11) 



The large behavior of the rj'g*g* vertex is the same 
as in the case of the gluon virtualities of the same signs 
[Eq. (O)] but involves a different function of the asym- 
metry parameter: 



ff2(w) hi^) ^-92 (-] = 



(4.12) 

The function g2{i^) is shown graphically i n Fig , o. In 
contrast to the function g2{'-^), given in Eq. ( I.IC ), which 
is antisymmetric in its argument, the function g2{'-^) is 
symmetric: g2{—uj) — g2{uj). 

When one of the gluons is on the mass shell, say, the 
second gluon (g| — 0), the leading gluonic contribution 
to the form factor is: 



F^%lg.(g^O,m^,) 



AnaM) Cql 



16B. 



(9) 



X < - 




.Ml 



(9) 



bB 



4(g2)2 



o^s{\ql\) 



In 1 



(4.13) 




In the limit of the large gluon virtuality (jg^j ^ m^,) the 
gluon contribution to the rj'g*g* vertex takes the form: 



4nas{\ql\) C 



16B. 



(9) 



k?l) 



bBk 



6 

(g) f as{\ql\) 



(4.14) 



IQl 
81 



As noted in the case of the quark contribution, when the 
gluon virtuality is time-like, q\ > 0, there is a threshold 



near qf — rn^, , but unlike the quark case, the gluonic 
contribution is regular. 



pia) 
v'g'g' 



(rn^,,0. 



Anasim^,) C 



16B 



(q) ( ojjr^^ 



5B. 



(4.15) 



as{^Jl) 



101 
81 



We point out that Eq. ( 4.13| ) is to be used for \qW> 



for the space-like gluon virtuality. In the region < 
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\qi \ < m^, the ry'-meson mass is the largest scale param- 
eter and in Eq. (4.13) the following change should be 
done: as((3^) ajjnl,). 



V. SUDAKOV-IMPROVED rj'g'g'' VERTEX 

In discussing the rj' g*g* vertex for the situation in 
which one of the gluons or both are far from the mass 
shell, the typical mass scale can be defined by the 
largest virtuality q\ or | ^ to^,). If the transverse 

momenta kj^.; of the partons in the ry'-meson are taken 
into account in the hard scattering approach, the per- 
turbative expansion of the vertex encounters large loga- 
rithms of the form ln(Q^/k^j), and it becomes manda- 
tory to sum the multiplc-gluon emissions. The formalism 
for the soft and coUinear gluon resummation was intro- 
duced by Collins and Soper |^ and Collins, Soper and 
Sterman ]3^ . Such gluon emissions give rise to pow- 
ers of double logarithms in each order of perturbation 
theory and their contribution exponentiates into the Su- 
dakov function ||3^ . The Sudakov exponents are known 
both for quarks from the Drell-Yan (DY) process and the 
deep inelastic scattering (DIS), and for gluons from the 
gluon fusion into 27, gauge, or Higgs bosons final states. 
This formalism is also suitable for the description of the 
hadronic wave-functions and the hadronic form factors, 
such as the electromagnetic and transition form factors 
of the pion ^ ^ ^ , and is also of interest here. 

In the modified Hard Scattering Approach 
(mHSA) the quark and gluonic invariant am- 

plitudes which define the rj' g*g* vertex can be written 
as: 



1 



Sp{75(f -m,Ori/)}, (5.1) 







(2^)2 



[T)!'U, (5.2) 



where b is the separation between the ?/ constituents 
in the transverse configuration space, often called the 
impact parameter, k_L is the transverse momentum of 
one of the constituents in the ry'-meson rest-frame, and 



the functions T^/^ and T]!^' are defined in Eqs. (|3.2| ) 
and (4.2), respectively. In this approach we take the 
Ty'-meson wave-function in the form similar to the pion 
wave- function 114, EQl : 



*(f)(x,Q,b) = 



27rC 



(cc, Q) exp 



xxh 



(5.3) 



where p = q for quark or antiquar k an d p = g for gluons, 
the constant C is defined in Eq. (2.3), and and 0^^^ 



have the forms presented in Eqs. (2.4) and (2.5), respec- 
tively. We assume that the transverse parts of the quark 
and gluonic components of the ry'-meson wave-function 
are universal and have a simple Gaussian distribution 
in the separation length b = |b|, i.e. ^^P^{x,Q,h) = 
^'■P^{x,Q,b). The parameter a can be determined from 
the average transverse momentum of the 77'-meson. For 
the numerical analysis we take the transverse size pa- 
rameter — 0.861 GeV, obtained for the 7r-meson by 
KroU et al. [|l^, |2^. The soft-gluon emission from the 
quark, antiquark and gluons inside the yy'-meson can be 
taken into account by inc luding the QCD Sudakov factors 
S'^P^{x,Q,b) in Eqs. and (Isll) . 

We shall perform our analysis by taking into account 
the next-to-leading logarithmic (NLL) contribution to 
the Sudakov factors, which are known for both the 
quark and gluonic ^ cases. In applications 
to form factors of mesons the Sudakov factor is ISStl : 



S<^P\x,Q,b) 



exp 



-Sp{xQ, b, Ci, C2) - Sp{xQ, b, Ci, C2) 



(5.4) 



where the Sudakov function Sp((5, 6, Ci, C2) for an external quark (j) = q) or gluon [p = g) line can be presented in 
the following general form p3[: 



Sp(Q,6, Ci,C2) = exp i -i j 



dq^ 



Aj,{as(/))\n^+Bp{asiq')) 
1 



(5.5) 



The integral term in Eq. (5.4) arises from the applica- tion of the renormalization group to the 77'-meson wave- 
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function as well as to the hard scattering amplitude. The 
upper limit is defined by the largest mass scale ap- 
pearing in the hard scattering amplitude. The functions 
7p(as) are the anomalous dimensions of the partons in- 
side the Ty'-meson in the axial gauge Efl : 



explicit expression for the LL and NLL terms of the Su- 
dakov functions in the form: 



SpiQ,b,C\,C2) 



A 



(1) 



00 



■Lc 



InA + i 
A 



(5.9) 



47r 



(5.6) 

Notice that as the energy parameter of the Sudakov 
function one should use the energies xQ and xQ of 
the ?7'-meson constituents which are, in general, differ- 
ent. The functions Ap(as) and B^ioLs) are perturbatively 
computable as power series in the strong coupling con- 
stant a,: 



4A, 



(2) 



B 



(1) 



/3g 



|-[lnA- A + l] + -^ In A 



PI 



(1) 



A In A + (1 - A) (1 + IuLq) -f Inig In A - - In^ A 



where A = Lg/Lb, Lq = ln{C^Q'^ / A'^), and Lb = 
ln(CiV62A2). 

Let us discuss the range of the impact parameter b. 
In general, it is defined as a positive variable but in 
the region b > 1/A the Sudakov exponent diverges due 

at 



— j Bp{as) — \~) p ' to Landau pole in the QCD coupling as{^P' 



(5.7) 

where A^'^ and are the expansion coefficients, spec- 
ified below for n — 1 and n — 2. We recall that the co- 
efficients A^p"^ are universal and lead to the resummation 
of the leading logarithmic (LL) contributions. The coef- 
ficients A^p'^ and Bp^^ give the NLL terms. These coeffi- 
cients are also process-independent but they depend on 
the renormalization and factorization schemes through 
the parameters Ci and C2. The next expansion coeffi- 
cients are process-dependent, as demonstrated recently 
for B^^^ [|l|. However, to the NLL accuracy, we do not 
require them. The coefficients of the LL and NLL terms 
of the Sudakov exponent for quarks Q and glu- 

ons M M are: 



^(1) 
q 

^(2) 

q 



^(2) 

g 



5(1) 



— Cp, 



X + /3o In 



2 



K 



-/3oln^' 



(5.8) 



l + 21n 



C2&0 
Ci 



6 In 



where Cp = (A^^- l)/(27Vc), Ca = {N^ = 3 in QCD), 
bo = 2exp(— 7^) where is the Euler constant, nj is 
the number of active quarks (with masses < 6*2(5), 
the coefficient K IBTI is: 



K = Ca 



67 
18 



2 1 



and Po has been specified earlier in Eq. (2^). Note that 



the coefficient Bq^^ given above in Eq. ( |5.8| ) is taken from 
Ref. jlj] and it differs from the corresponding coefficient 
given in Ref. 

The knowledge of the LL and NLL perturbative coef- 



ficients A^p \ and Bp^' as well as the strong cou- 

pling as{q^) at the two- loop level (^) allows to get the 



?(i) 



ji = A. In this case the perturbative calculation is no 
longer valid and a prescription for parametrizing the non- 
perturbative physics in the low transverse momentum (or 
large b) region is necessary. We shall restrict ourselves to 
the region b < 1/A. There also exists a lower bound 
on the Sudakov function, arising from the consideration 
that the transverse momentum of the emitted gluons is 
not allowed to be large, i.e., we have to restrict to the 
region |k^| ^ C2Q, where C2Q is the hard scale. Thus, 
the region of applicability of the Sudakov formalism is 
limited to the region: I/C2Q < 6 < 1/A. 

In the following it is necessary to fix the arbitrary pa- 
rameters Ci and C2. The constant Ci determines the 
onset of the non-perturbative physics. The renormaliza- 
tion constant C2 specifies the scale of the hard scattering 
process. Following Ref. [^, we make the choice Ci — 1 
and C2 = I/V2 in our analysis ||4^. 

Using the definitions of t he q q and gluonic parts of the 



r?'-meson vertex ( |3.4| ) and ( [4.4[ ), respectively, the results 



for the corresponding vertices, including the dependence 
on the transverse momenta, are: 



i^^'j.,.('zf,92,K') = 4™,(g^)^ 
1 

X JdxJ^¥'''>{x,Q,h) S^''\x,Q, 



(5.10) 



b) 



{2n) 



-L g-z(k_Lb) 



: \- {x ^ x) 



xqf + xq2 — xxm^, — + 



dxj'^^¥^\x,Q,h)S'-'\x,QM 



(5.11) 



-»(k_Lb) 



ml, + {x-x){ql-ql) 
xq\ + xq2 — xxm^, — -|- ze 



(a; <-> x) 



For the space-like gluon virtualities [qf < 0, g| < 0, 
and = —{qf + (j|)] the integration over the transverse 
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momentum kj^ can be done resulting in the modified 
Bessel function of order zero, Kq{z). Taking this int o ac - 
count and substituting the ?7'-meson wave-function (^.3[) , 



J 



the quark and gluonic contributions reduce to the two- 
dimensional integrals of the form: 



1 oo 

F^?l,g,iQ,u;,7^) = -47:a,{Q')-^J dx^<-''\x,Q) j dhh^-^^''" 1^'^" S'^^\x,Q,b) Ki+\x,hQ) 



(5.12) 



= -4™,(g2) 



C 



dx (j)'''^\x,Q) I c?6a6a exp 



'4a2A2 



i^,;?].g.(Q,c^,77) = -4™,(Q2)c ( dx<jy^a\x,Q) /" d& fee"^^"'/*"' Q, 6) 



hi 14-^ {x, bQ) - [x ~x)u: K^+^ (x, bQ) 



(+), 



-4^a,(Q2)^ I dx^^<^\x,Q) I dbt^bt, exp 



4a2A2 



"A 



X S^3\x,Q,b) \hK^^ \xMQa)-{x-x)uoK^^\xMQa) 



where 



S'^''\x,Q,b)K^+\xMQA). 

(5.13) 



Kl,^\x,bQ)^K^^\xMQK) = -[KoibAQAX+ix,u;,v)) ± KoibAQAX+{x,iu,T]))] 

I 



with the dimensionless parameters defined as 6a — bA, 
Qa = Q/A, \r]\ = the function A|(x,cj,7y) = 

[1 + uj{x — x) ± 2xx\'r]\]/2, and the quark and gluonic 
components of the ry'-meson wave-function are taken in 
the form of Eq. (|5.3|). The restriction 6 < 1/A due to the 
Sudakov formalism is also taken into account in the final 
result for the space- like gluons presented in Eqs. ( 5.12| ) 
and ( |5l3| ). 

The arguments of the Bessel functions 
KQ(bQX+{x,LL!,r])) and Ko{bQX+{x,LL!,r])) in the r/'- 
meson vertex de fine the largest scale parameter t in the 
Sudakov factor (5.4): 



t — max {QX^{x, uj, rj), Q\^{x, uj, rj), 1/6} . (5.14) 

The dependence of the Sudakov factors S^'^^{x,Q,b) 
and S*'^' (a;, Q, 6) on the impact parameter b at x — 1/2 
and Q = lOA, 20A and 50A is presented in Fig. ^. It is 
seen that the Sudakov factors give a cut-off in addition 
to the restriction 6 < 1/A of the cut-off of the integration 
interval over 6. With increasing the typical mass scale Q 
this effect is stronger in the quark contribution. 

In Fig. ^ we show the distribution of the quark and glu- 
onic parts of the vertex in the impact parameter space 
integrating here with a variable cut off, bcut- The curves, 
showing the dependence of qf F^fg,g{ql,0,m'^,) on bcut, 
start from zero for bcut = and reach their full height and time-like regions {as{—q'^) = as{q'^)[l+0{as)])- The 



for bcut — 1 beyond which we consider any remaining 
contributions as non-perturbative. The faster rise of the 
rj' g* g* vertex is observed for the quark contribution as Q 
increases while the gluonic one is mildly sensitive to the 
typical mass scale. At the largest value shown, Q = 50A, 
the quark curve is quite flat for bcut > 0.2 which indi- 
cates a small contribution from this region. As for the 
smaller values, for example, Q — lOA, the vertex gets 
a rather large contribution from large distances in the 
impact parameter space. The gluonic part of the vertex 
receives the main contribution from the short distances 
6 < 0.3. In contrast to the quark contribution, the glu- 
onic contribution is negative and decreases the summed 
vertex function by approximately 20%. 

To obtain the jy'-meson vertex function for the time- 
like gluon virtual ities one h as to replace (32 _q2 
\ri\ —rj in Eqs. ( 5.12 ) and ( 5.13 ). The vertex function is 
an analytic function of the scale parameter and it can 
be analytically continued from the space- like region to 
the time-like one. As mentioned earlier, the rj' g*g* vertex 
results from convoluting the 77'-meson wave-function and 
the hard scattering amplitude kernel. As we do not take 
the QCD corrections to the hard scattering amplitude 
into account, the strong coupling, as{q^), has the same 
dependence on the scale parameter in both the space-like 
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FIG. 8: The Sudakov function plotted in terms of the di- 
mensionless impact parameter 6a for quarks (upper figure) 
and gluons (lower figure) with a; = 1/2 and Q = lOA, 20A, 
and 50A. 



FIG. 9: The quark-antiquark (upper figure) and gluonic 
(lower figure) contributions to the r]'g*g vertex function 
(p) (r3 n versus the cut off in the impact parame- 



ter space, bcut, for the space-like gluon virtuality 
with Q = lOA, 20A, and 50A. 



other component of the ?7'-nieson wave-function sensitive 
to the transition from the sp ace- like region to the time- 
like one is the Sudakov factor ( ^.4| ). The analytic continu- 
ation of the Sudakov factor was analyzed in Refs. [^6|, Q 
with the result that the time-like and space-like form fac- 
tors have the same scale dependence in the asymptotic 
regime. We use the space-like expression for the Sudakov 
factor to study the rj' g*g* vertex for the time-like gluon 
virtualities. The hard scattering amplitude kernel gives 
the additional transfer momentum dependence due to the 
quark and gluon propagators. We recall that, in gen- 
eral, the propagator is defined in the complex momentum 
plane and, hence, it can have both real and imaginary 
parts in the time-like region of the momentum due to ie 
term in the denominator. Here, we perform the analy- 
sis of the propagators in the transverse momentum space 
keeping their transverse momenta dependence. In mHSA 



it is more convenient to operate with the Fourier trans- 
form of the propagators in the impact parameter space 
but this analysis can be easily reformulated in terms of 
the impact parameter b, which is the conjugate variable 
to kj_- 

The integration plane can be divided in several re- 
gions as shown in Fig. |l^. The interval < |kx| < A is 
the non-perturbative region. We note that the Brodsky- 
Lepage approach, in which the dependence on the trans- 
verse momentum is neglected, corresponds to the point 
kj^ = in Fig. |l^. It is seen that this transverse mo- 
mentum value lies in the deeply non-perturbative region 
and too far from the contributing region of the transverse 
momentum integration. It means that in the Brodsky- 
Lepage approach the imaginary part of the r]'g*g* vertex 
can not be calculated correctly and the natural prescrip- 
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FIG. 10: Different regions in the integration plane in 
the variable k^. The abbreviations used correspond to 
non-perturbative (NPR), soft-gluon (SGR), and hard-gluon 
(HGR) regions. The point |kxj = corresponds to the 
Brodsky-Lepage approach. 



tion is to drop it altogether. Another consequence of the 
approximate forms of the propagators in the Brodsky- 
Lepage picture is the appearance of the singularity in the 
r]'g*g* vertex in the region close to the Ty'-meson mass. 
The interval jk^j > A can be effectively divided into two 
regions: the soft-gluon region (SGR) where the Sudakov 
factor is taken into account and the hard-gluon region 
(HGR) where the Sudakov effect is absent. In Fig. |l^ 



we also point out the value of v^ia" to show the im- 
portance of the transverse momentum distribution of the 
r^'-meson wave-function. In the region |kj_| > ^/xEa~^ 
the contribution to the vertex is effectively suppressed 
by the assumed Gauss ian behavior of the transverse mo- 
mentum distribution (5_^). 

To calculate the integrals (5.1C) and ( p.ll ) over the 
transverse momentum k_L in the time-like region of the 
gluon virtualities it is useful to do the analytic contin- 
uation of |k_L| in the complex plane. The new feature 
with respect to the space-like case is that the contour of 
the transverse momentum integration goes now near the 
pole located at k^ = + ie. This pole, except in 



the end point regions [a; ^ for the second gluon on the 
mass shell, or a; ^ 1 for the first one due to the func- 
tion A_ (x, w, 77)] , is far from the bounds of integration of 
the variable |kj^|. Therefore, this integral can be evalu- 
ated by deforming the contour of the integration in the 
complex plane of this variable. 

The pole encountered in our computations when an 
internal quark or gluon goes on the mass shell do not 
correspond to an observable state. They appear in the 
effective t]' g* g* vertex which is not an observable by it- 
self. Provided this vertex is included in the complete 
amplitude of some real process such as -B — > rj' K and 
B rj'Xs, this amplitude will have an additional com- 
plex phase factor as a reminder of this singularity. Notice 
that the complete physical amplitude of the process in- 
cluding the effective ii g* g* vertex in a purely hadronic 
computations has also the poles reflecting the existence 
of the intermediate physical (mass-shell) states which are 
hadronic ones. Therefore, it is natural to expect the ap- 
pearance of the additional phase factor here also. 

According to the arguments given above, the rj' g* g* 
vertex for t he tim e-like gluo n virtualities can be obtained 
from Eqs. ( 5.12 ) and ( 5.13| ) by changing \rj\ —rj and 
A'o(6(5A+(a;, w, 77)) Ko{ibQX-{x,uj,ri)). If one uses 



the relation Ko(iz) 



-iiT H, 



(2) 



(z)/2 where 



(2) 



{■A 



Jo{z) — ilo(z), (here, H^'^\z), Jo{z) and Yq{z) are the 
second Hankel function, Bessel function of the first kind, 
and Neumann function, respectively |Q), it is seen that 
in contrast to the case of the space-like gluon virtualities 
the vertex gets an imaginary part. After these changes 
the quark and gluonic contributions to the rj' g* g* vertex 
function are: 



J 



exp 



'4a2A2 ^. 



2A2 



dx4)'^3\x,Q) / dbAbA 



exp 



'4a2A2 



S^''\x,Q,b)H^+\x,bAQA), (5.15) 

(5.16) 



X 5(s)(x,Q,6) Ij^H^o \x,bAQA) + ix-x)ujH^o'^\x,bAQA) 



where 



H^o^\x,bQ) = H^'^\x,bAQA) = l \h^^^ {bAQAX-{x,u;,v))±H^o'^ (6aQaA_(x, c^, 77)) 



The real and imaginary parts of the rj' g*g* vertex in the 
form q\ J^^fg.g are shown as functions of the impact pa- 
rameter cut off, 6cut> in Fig. |ri|. As in the case of the 
space-like virtualities, a fast rise is observed with Q in 
the region of small values of the impact parameter cut-off. 



The oscillatory behavior of the cylindric functions results 
in an oscillatory form of the vertex function for small and 
moderate values of the impact parameter. The contribu- 
tion to the vertex from the large distances, b > 0.8, is 
small. It is also seen that as the scale parameter Q in- 
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FIG. 11: The real (solid curves in the upper figure) and imag- 
inary parts (dashed curves in the upper figure) and the ab- 
solute value of the vertex function qiF^i g* g{q\,Q, rn^,) (lower 
figure), plotted as functions of the cut off, bcut, in the impact 
parameter space for the time- like gluon virtuality, qi ~ Q^, 
with Q = lOA, 20A, and 50A. 



creases, the relative contribution of the imaginary part 
in the absolute value of the vertex is strongly suppressed, 
and the vertex is mainly defined by the real contribution 
in the large asymptotics. 

The asymptotics of the quark part of the vertex F^f'g-g^ 
in the limit of large is: 



- 4..(g^) 



3C 



f^{u) + — y+y_ (5.17) 



X (1-y+y-) ilf (QAj;+)-i?r(QAy-) 



(2), 



where y± = ^(1 ± cj)/2, and the function f{){uj) is de- 
fined in Eq. (3.11). In this expression the large argu- 
ment asymptotics of the Hankel function should be used: 
H^^\z) ~ exp(-^[z - 7r/4]) @. The same cor- 

rection of order to the gluonic part of the vertex 

function Pj^fg^g-^ is equal to zero and one can use the 



asymptotic behavior defined by Eq. (4.9). 



VI. NUMERICAL ANALYSIS AND 
COMPARISON WITH EXISTING RESULTS 

In this section we give a numerical analysis of the 
vertex in the time-like and space-like regions, de- 
rived in the preceding sections. To that end we spec- 
ify the input parameters. The dimensional factor C 
[see Eq. (^)] contains the decay constants fq and fs 
and the mixing angle which can be constrained from 
the existing experimental data. We shall adopt here 
the Feldmann-KroU-Stech rj — rj' mixing scheme [ pO| , 
which is phenomenologically consistent and satisfies the 
constraints from chiral perturbation theory [ p^ , yield- 
ing @: /g = (1.07 ±0.02)/,, /, = (1.34 ± 0.06) A, and 
(j> = 39.3° ± 1.0°, where ~ 131 MeV is the pion decay 
constant. Using the central values of the parameters, one 
gets C ~ 2/, ~ 260 MeV. 

The quark and gluonic wave-functions (Al) and ( [A^ ) 

contain both free (i?!''' , Bn ^ ) and constrained (pi''' , pri ) 
parameters. The constrained parameters depend on the 
anomalous dimensions and are given in Appendix 
while the free parameters can be fitted from the ex- 
perimental data, for example, from the r?'7*7 transition 
form factor. We take the following restrictions on the 
first correction to the leading order wave- function (|Al|): 

IB^"^^] < 0.1 and \pi^^ B!f^\ < 0.1 in order to keep 
(l)^'^\x, Q) close to its asymptotic value: (f>a.s{x) = 6xx, in 
agreement with the experimental data on the rj'j*^ tran- 



sition form factor |15[. Taking into account P2^^ ~ —1/90 
from Eq. (|A10|) we get||| {B^^^'^l < 9.0. Below, we shall 



present the rj' g* g* vertex function for the maximum al- 
lowed values of the non-perturbative parameters, i.e.. 



\B 



(9) I 



0.1 and \B. 



(9) I 



9.0. We have also studied nu- 



merically the resulting vertex function for considerably 
smaller values of these parameters. In that context we 
note that as the NLO quark contribution to the overall 
T]'g*g* vertex function is small, there is not much sen- 
sitivity to the variation of the parameter i?2'''- Hence, 
we fix this parameter to its maximum allowed value 
|i?2'''| = 0.1. However, there is considerable sensitivity to 
the variation of the parameter Bj^'' . To show this we shall 
take 1-82"'' I — 3.0. The resulting theoretical dispersion be- 



tween the two cases 
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FIG. 12: The rj'g'g vertex J;,/^ 



of with Bj'' =0.1 and two values of B, 



{qi , 0, m^/ ) as a function 



5(9) 
-'2 



= 9.0 (the 



upper plot) and Bj^'' = 3.0 (the lower plot) in the Brodsky- 
Lepage approach. The dashed curves are the leading (LQ), 
next-to-leading quark-antiquark (NLQ), and gluonic (NLG) 
components, and the solid curve is the sum. 



be taken as representative of the theoretical uncertainties 
due to the undetermined non-perturbative parameters. 

First, we consider the case when one of the gluons is on 
the mass shell, which we take for the sake of definiteness 
to be the second one (gl = 0). In Fig. ^ we show the 
leading (long-dashed curve) and next-to-leading (middle- 
dashed curve) quark contributions as well as the gluonic 
one (short-dashed curve) for the time-like gluon virtual- 
ity (qf > m^i) corresponding to the maximal values of 



the free parameters: b'^'^ =0.1 and B^' — 9.0 (upper 
plot), and for B^''^ = 0.1 and b[^'^ = 3.0 (lower plot). The 
solid curve in each of these figures is the total contribu- 
tion to F^'g.g- (g^, 0, m^/) in next-to-leading order. One 
can see from this figure that with the given parameters, 
the next-to-leading quark contribution is subdominant 
in the entire region of shown, except in the neighbor- 
hood of the threshold ~ m^y , where both the leading 



FIG. 13: The r)' g*g vertex for the time-like gluon virtual- 
ity with B. 



0.1 and two values of B. 



(9). 



B. 



(9) 



9.0 



(upper solid curve) and Bj^' — 3.0 (lower 
calculated in the Brodsky-Lepage approach in 
The long- and short-dashed curves correspond 
tions F^fp^{qlO,ml,) = 1.8 GeV-V(g?/m^, 

F^"PgiQ^) = VSas{Q^)/i7r^) i, respectively. 



form of the vertex function F^f'V — 1.8 GeV 

§'7'9*9 
] is also shown as a dashed line. 



solid curve) 
this paper, 
to the func- 
- 1) i and 
A constant 
suggested in 



and non-leading quark contributions have logarithmic di- 
vergences. Notice that in this region our results are not 
valid and we can trust only the region 3> m^, . The 
gluonic contribution (called NLG) for the maximum val- 
ues of the free parameters is comparable to the leading 
quark contribution (called LQ), as shown in the upper 
plot in Fig. |l^, increasing the rj' g* g vertex by almost a 
factor 2 as compared to the case when only the quark 
content of the 77'-meson is assumed. This may be con- 
sidered as the maximum gluonic content of the ry'-meson 
allowed by current data. Even in the more realistic case 
with bI^^ =0.1 and B^'' = 3.0, we see from the lower 
plot in Fig. [ij that the gluonic contribution is not small, 
and also in this case it enhances the value of the total 
9*9 vertex function (the solid curve in Fig. |l^) at the 
level of few tens percent. 

Let us compare the if g* g vertex calculated here with 



(<?) _ 



the two set of values B.^ 
(the upper solid curve in Fig. 



0.1 and B. 



(a) 



B. 



(9) 



-"2 

and S^"*' - 



9.0 



3.0 (the lower solid curve in Fie 



= 0.1 and 
jj) with the 

ones given in the literature. In Ref. M, Kagan and 
Petrov have parametrized the if g*g vertex in the form: 
F,yg.g{qlQ,ml,) = HiO,0,m^^,)/{qf/m^^,-l), where 
H{0,0,m^,) ~ 1.8 GeV~^ is a phenomenological param- 
eter extracted from the experimental data. Its behavior 
is described by the long-dashed curve in Fig. |l^, labeled 
as F^^^gfg . In Ref. |^ , Hou and Tseng have parametrized 
the ?]'g*g vertex as: Fjj>g.g{qf) = V3as{qf)/{nf^), which 
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FIG. 14: The r)' g* g vertex function qiF^'g*g{qi) for the 
space-hke gluon virtuahty with Bj'^ = —0.1, Bj®' = —9.0 
and Bj''' ~ 0.1, Bj^' = 9.0 (the upper and lower solid curves) 
and B;^"' = -0.1, B^''' = -3.0 and B^"' = 0.1, B^"' = 3.0 
(the upper and lower dotted curves) in the Brodsky-Lepage 
approach. The dashed curve is the leading order contribution 
to the function obtained by setting Bj'' = Bj"' — 0. 

is drawn as the short-dashed curve in Fig. O, labeled as 
^v'g'g' -f'i'^^lly' ^l^*^ show the constant vertex func- 
tion Fn'g-.g{ql) = F^fg^J^ 1.8 GeV"\ assumed by At- 
wood and Soni in Ref. H. It is seen that our result (the 
upper solid curve) compares well with the form used in 
Ref. for a wide range of q\, but is significantly differ- 
ent (smaller) than the one assumed in Ref. ||] and is in 
complete disagreement with the one assumed in Ref. . 
Our analysis shows that in the case of the lower pos- 
sible values of the free parameters: B^2^ = —0.1 and 
b'^^ = —9.0, the vertex function Fri'g'g*{qf,0,m'^,) de- 
creases and a substantial disagreement with the vertex 
suggested in Ref. would result. However, as discussed 
below, data on the electromagnetic transition form factor 
of the yy'-meson disfavors this choice of the parameters. 

We show the results of our calculations for the function 
qf Fri' g" g((i\) in the space-like region of the gluon virtual- 
ity = -Q^ < 0) in Fig. |l|. The dashed curve in the 
middle is the leading contribution to the rj' g*g vertex, 
obtained by setting the non- asymptotic parameters 
and B^'^ equal to zero. The upper and lower solid curves 
are the vertex functions F,yg*g{qi, 0, m^,) taking into ac- 
count the next-to- leading correction with B^'^ = —0.1, 
b''^^^ = -9.0 (upper curve) and B^''^ =0.1, S^^^ = 9.0 



(lower curve). The dotted curves correspond to the choice 
= -0.1, B^^^ = -3.0 (upper curve) and b'^^ = 0.1, 
B^2^ =3.0 (lower curve). We note that the correction is 
mainly determined by the value of the parameter B^'' . 
It is seen that at small values of the gluon virtuality 
{Q^ < 5 GeV^) the complete correction becomes approx- 
imately equal to or even larger than the leading order 
contribution. This implies that in the region < few 
GeV^ our approximation is not valid and the contribu- 
tions of the next higher order corrections become impor- 
tant. The behavior of F,,'g.g(Q^, 0, to^,) for large is 
in qualitative agreement with the corresponding electro- 
magnetic transition form factor, F^/^.^(Q^, 0, m^,), mea- 
sured in 77*-collisions, more recently by the CLEO and 
L3 collaborations |l5). Hence, this would suggest that 
for the free parameters we should use the positive val- 
ues: i?2''' = 0.1 and 52^'' = 3.0, or values close to them, 
yielding the lower dotted curve in Fig. |l^. 

We shall now present the effect of the transverse 
momenta on the vertex functions. To that end, 
the ri'g*g* vertex functions \F^rg*g,{qf,q2,m'^,)\ and 

9^ -Fj)'g*£(* (<?! , ™^') presented in the Brodsky- 

Lepage approach in Fig. |l5|, for the space-like gluon vir- 
tualities with q^ = qf + g| and no transverse momen- 
tum effects taken into account. The corresponding vertex 
functions in the mHSA formalism are shown in Fig. 
in which transverse momentum effects are included as 
discussed earlier. The various curves shown correspond 
to the following values of the second gluon virtuality: 
= 0, -1, -5, -10, and -25 GeV^. We have fixed 

-82'''' — 0.1 and -82^'' = 3.0 from the analysis of the Fr/'-y-j 
form factor. Both approaches give the same asymptotic 
behavior at large values of the mass scale of the form: 
\F,fg.g* \ ~ 0.4 GeV/Q^. In the region of small Q^, the 
mHSA-prescription modifies the behavior of the vertex 
functions as it decreases the gluonic contribution. 

The results for the r]'g*g* vertex with the time- 
like gluon virtualities {qf > 0, > 0) in the 
forms F,,/g-g.(g?,(7|,TO2,) and q^ F,y g. g. {qf , q^ , m'^,) with 

B;^"^' = 0.1 and B^f'^ = 3.0 in the Brodsky-Lepage ap- 
proach and mHSA are presented in Figs. ^ and |l8|, re- 
spectively. We recall that in the considered region of the 
virtualities, an imaginary part appears in the rj' g* g* ver- 
tex in the mHSA formalism in contrast to the Brodsky- 
Lepage scheme in which the rj' g* g* vertex is defined as 
real by prescription (see Sec. In both approaches, 
the absolute values of the function describing the r]'g*g* 
vertex decrease with increase of the second gluon virtu- 
ality, gf. When the virtualities of both the gluons are 
close to the Ty'-meson mass squared an enhancement of 
the vertex function is observed in both approaches. Tak- 
ing into account the transverse momentum distribution 
of partons inside the ?7'-meson in the mHSA, the singular- 
ity in the form factor obtained in the Brodsky-Lepage ap- 
proach is replaced by the resonance-like behavior. Never- 
theless, qualitatively the behavior of the vertex function 
in the Brodsky-Lepage approach (upper plot in Fig. [T^ ) 
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FIG. 15: The rj'g*g* vertex functions for the space-hke gluon 
virtualities in the forms: Ifij/g.g* (gf , q|, m^/)| (upper figure) 



and q^Fjfig-,g,{qi,q2,m^^,) (lower figure) at Bj'^ = 0.1 and and q^ Frjig- g* {qi,qi,mP) (lower figure) with Bj''' ~ 0.1 and 



FIG. 16: The r)' g* g* vertex functions for the space-like gluon 
virtualities in the forms l-Fjj'g^g* (qi, 92 , )| (upper figure) 



(9) 



B':^^ = 3.0, calculated in the Brodsky-Lepage approach. The 3)^' = 3.0, calculated in the mHSA formalism. Legends are 



(9) 



legends are as follows: = (solid curve), = —1 GeV^ 
(long-dashed curve), = — 5 GeV^ (medium-dashed curve), 
-10 GeV^ (short-dashed curve), ql = -25 GeV^ (dotted 
curve) . 



ql 



the same as in Fig. 15 



VII. INTERPOLATING FORMULAE FOR THE 
n'g'g* VERTEX 



and in the mHSA (left-down plot in Fig. ^8|) is similar. 

The typical values of the gluon virtualities at which 
the vertex function is well described by its asymp- 
totic behavior F^igtg* ~ 1/Q^ can be determined from 
Fig. |l^ in both the Brodsky-Lepage and mHSA ap- 
proaches for the case when one of the gluons is on 
the mass shell. It is seen that for large the ver- 
tex functions in both of these approaches are well cor- 
related. In the case of the negative gluon virtuali- 
ties, the function i^r;'g*g* reaches its asymptotic form 
at smaller values of (Q^ ~ 5 — 10 GeV^) than in 
the case of positive virtualities [Q^ ^ 50 — 100 GeV^). 



For the applications in various decay and production 
processes it is useful to find some approximate formulae 
for the vertex functions which are simple and can be used 
over a large domain of the gluon virtualities. We recall 
that Brodsky and Lepage |Q presented an approximate 
form for the tt — 7 transition form factor which interpo- 
lates between the PCAC value and the QCD prediction 



in the large Q region. Subsequently, in Ref. 



this 



form was extended to the case of the rj — ^ transition 
form factor. Very much along the same lines, a similar 
expression can be written for the ii g* g* vertex: 



2V3UD{q\u;) 



(7.1) 
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we use 



The above expression reproduces both 



the anomaly value and the large asymptotics of the 
vertex function: 



2 ■ -^1 - Ana (n?,\ ^ 

Q2^Q - 47r2J^' 



(7.3) 
(7.4) 



.=> 0.25 



0.00^ 



Presented in the form given in Eq. (7.1), the vertex 
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FIG. 17: The r)' g* g* vertex functions F^>g'-g*(qi,<lh^^') 
and q^Fnigtg-, {qf, q^, rn^,) for time-like gluon virtualities with 



F^]^, {q^ , oj) is a smooth function in the space-like re- 
gion of the gluon virtualities but has a pole at — 
Stt^ f'^D{q^ ,u}) I y/i in the time-like region. A similar be- 
havior for the Ti'g*g* vertex in the time-like region was 
obtained by Kagan and Petrov |Q from the evaluation of 
the triangle diagram. 

The mHSA approach removes the unphysical singular- 
ity from the vertex in the time-like region of the gluon 
virtualities but the vertex gets an imaginary part. In 
this case the real part of the approximate formula inter- 
polates the anomaly value and the large asymptotics, 
while the imaginary part goes to zero as — > 0. It is 
possible to get the approximate formula F^^?^{Q'^ ,uj) 
without the pole behavior in the time-like region of the 
gluon virtualities, if one takes into account the correc- 
tions of order 1/Q^/^ ( ^.17 ) to the vertex f unct ion. In 
this case the functional form defined in Eq. (7.1) is still 
valid but with the replacement: 



2ni 



■y+y- 



(7.5) 



B^''' = 0.1 and B^^' = 3.0 in the Brodsky-Lepage approach. tion H^f' {Q/^y±) is assumed in the limit of the large 



,(9) 



X (l-y+y-) H^^\QAy+) - H^^HQav-) 



where y± is defined by Eq. ( ^.17| ) and the Hankcl func- 

(2) 



The legends are as follows: (?2 = (solid curve), ^2 = 1 GeV 
(long-dashed curve), g2 = 5 GeV^ (medium-dashed curve), 
g| = 10 GeV^ (short-dashed curve), g| = 25 GeV^ (dotted 
curve) . 

where the largest energy scale parameter — q^ for 
the time-like total gluon virtuality and = —q^ for the 
space-like one. We introduce also the following function: 



asymptotics of it argument. As the function 'D{Q'^,uj) 
is complex, the vertex F™^^^(Q^, uj) becomes a complex 
function as well. In t he limit of large , Eq. ( [7.l|) with 
the replacement (7.5) reproduces the asympto tics (5.17) 
of the vertex function. The real part o f E q. (7.1) gives 
the value determined by the anomaly (7.3) in the limit 
—> 0. As for the imaginary part of the vertex func- 
tion, it goes to zero as ImF^/^^^(Q^, w) ^ Q"* in the 
limit of small virtualities according to the approximate 



D{q^,u;) 



+ gl.92(w) 



(7.2) 



formula (7.1 



(9) 



5B}, 



(s) f as{Q^) 
as(Aio) 



When one of the gluons is on the mass shell (g| = 0) the 
asymptotic functions /o(^) and g2{'-^) have the following 
values: /o(l) = 1 and 32(1) — 1/6. We do not include 
the next-to-leading quark contribution which is strongly 
suppressed. In the numerical analysis the usual prescrip- 
tion for the evaluation of strong coupling should be used: 
it is estimated using the value of the largest scale param- 
eter of the problem. For the anomaly, with real gluons. 



VIII. SUMMARY 

In the present paper we have analyzed the rj' g*g* ver- 
tex, Fjj' g* g* {qf , qf , m'^,) , for off-mass-shell gluons in the 
hard scattering approach. The ?7'-meson wave-function is 
evolved using the evolution equations for the quark and 
gluonic components. In difference to Ref. it is shown 
that, within the possible variation of the parameters 62'^^ 

and 82^^ of the r^'-meson wave-function, the gluonic con- 
tribution can not be ignored. To further quantify this. 
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FIG. 18: The ri'g*g* vertex F^j' g* g* {qf , qi , m^,) in the mHSA formahsm, with the time-hke gluon virtuahties, — 0-1 s-nd 
52"' ~ 3.0 where = 9i + 92- Legends are the same as in Fig. 



one needs to know the non-perturbative parameters en- 
tering in the evolution of the r/'-nieson wave-function. Us- 
ing the maximum values of these parameters |i?2'''' | = 0.1 

and l^^^^l = 9.0, we find that the gluonic contribution is 
comparable to the quark contribution, leading to a siz- 
able (almost a factor 2) enhancement of the rj' g*g* ver- 
tex. The magnitude of the vertex is not very sensitive to 
the value of the parameter b'^^ , and even a factor three 
to five reduction from its default value used here does 
not change the vertex function appreciably. However, 
the vertex function is rather sensitive to the input value 
of the parameter B':f^ . We have explicitly shown this de- 
pendence in our work. Hence, it is important to constrain 
or measure this quantity for more definitive predictions. 
However, we find that even with a significantly reduced 



value of say for |i?2^''| = 3.0, the gluonic correction 
is of the order of a few tens percent, and hence it must 
be included in all perturbative treatments of the if g*g* 
vertex function. 

We have obtained analytic expressions for the vertex 
describing the ry'-meson transition into two gluons with 
arbitrary virtuahties in the Brodsky-Lepage approach in 
which the transverse momentum dependence of the par- 
tons inside the 77'-meson is ignored. The appearance of 
a singularity in the region of the Ty'-meson mass indi- 
cates that this approach is valid only in the asymptotic 
region for the total gluon virtuality, = (?i + q|, i.e., 
for \q^\ ^ 771^,. In this region the r]'g*g* vertex has 
the usual behavior as the pseudoscalar meson form fac- 
tors: Fj^'g'g' ~ l/?^- We have compared our results 
with some of the existing parametrizations of the rj' g*g* 
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19: The large asymptotics ol the vertex functions 
g(gi,0, m^,) with the time-like (+) and space-like (— ) 



gluon virtualities with B. 



0.1 and B. 



(9) 



3.0 in the 



Brodsky-Lepage (BL) and the mHSA approaches. 



vertex, Frfigrg* {q\, 0, m^,), in the time-like region. Corre- 
sponding results for the vertex F^/g.g. (g^, 0, m^,) in the 
space-like region of are also presented, and the be- 
havior of the T]' g*g* vertex is found to be close to that 
of the electromagnetic transition form factor of the ry'- 
meson - an information which has been used to fix the 
parameters B^^ and B^'' in the evolution of the ry'-meson 
wave- function. 

The mentioned singularity in the Brodsky-Lepage ap- 
proach can be circumvented in the modified hard scat- 
tering approach (mHSA), in which the transverse mo- 
mentum dependence of the hard scattering amplitude as 



well as the transverse momentum distribution and soft- 
gluon emission (the Sudakov effect) in the yy'-meson wave- 
function are taken into account. Analytic properties of 
the rj' g* g* vertex are reviewed to obtain the correct ex- 
pression in the time-like region of the gluon virtualities. 
It is shown that due to the ie prescription of the prop- 
agators, the rj' g* g* vertex acquires an additional phase 
factor in comparison with the space-like expression. We 
present the resulting rj' g*g* vertex in the mHSA for- 
malism for the space-like and time-like virtualities and 
work out the large asymptotics. Numerical analysis 
shows that in the modified hard scattering approach the 
rj' g* g* vertex in the time-like region of the gluon virtual- 
ities reaches the asymptotic form at rather large values, 
Q > 10 GeV, while in the space-like region, the corre- 
sponding vertex can be used in its asymptotic form al- 
ready at Q ^ 3 — 5 GeV. Finally, we have provided simple 
interpolating formulae for the vertex F,y g<- g<. (q'^ , to) for 
the time-like and space-like gluon virtualities, which re- 
produce the anomaly (for on-shell gluons) and the asymp- 
totic form, determined in the hard scattering approach. 
The computed vertex F^ji g* g* {Q^ , lo) has ready applica- 
tions in a large number of decays and production pro- 
cesses involving an ry'-meson. 
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APPENDIX A: SOLUTIONS OF EVOLUTION 
EQUATIONS 

The solutions for the quark and gluonic wave- functions 
of the color- and SU(3) p flavor-singlet pseudoscalar me- 
son were obtained in Ref. Their general 
forms are: 



{x, Q)=&xx\l + ^ 



X 



(Al) 
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<l)^s\x,Q) =xx J2 



even 7i>2 



3(9)5(9) 



(A2) 



Here, is the scale of the hard process, /xq — 0.5 GeV 
is the typical hadronic energy scale below which no per- 
turbative evolution takes place, x and x = 1 — x are the 
momentum fractions of partons inside the pseudoscalar 
meson, (x) are the Gegenbauer polynomials of the or- 
der n with the index v |30[ . When the index is semi- 
integer, h' — m + 1/2, the Gegenbauer polynomials can 
be defined by the following relation: 



2'"to! 
2"n!(2m)! V^^^ 



m-\-n 



{x^ ~ 1)". (A3) 



In our analysis we are interested in the polynomials of 
the first three orders: 



-1/2 



(x) 



1, 



(2m -I- \)x, 



(A4) 



C2"+'/'(x) i(2m+l)[(2TO + 3)a;2-l] 



Hence, C'^'^[x) 
differential equation: 

d 

dx 



|(5a;2 



1) and Cl^'^ix) = 5a;. The 
C:{x)^2,yC:tl{x), (A5) 



allows to connect the coefficients of the highest powers of 
these polynomials — {3/n)^/^^an-i, where the 

following representation for the polynomials is assumed: 

C^x) = x" + >„^i -I- . . . + >o- (A6) 



Eqs. (Al) and (A2) contain the set of parameters 
called 75, defined as: 



2 



where the anomalous dimensions 7,^- are ||2^, |2J, |2£ 



.(A7) 



7gq 
7qg 



Cf 

Po 

Cf 

Po 3(n+ l)(7i + 2)' 
Uf 12 

(n + l)(n + 2)' 



n-l-l 



(n + l)(n + 2) 
n(n + 3) 



^^7 



(A8) 



7gg 



a 1 
? 4Vi 

(n-|-l)(n + 2) 



Here, Nc is the number of colors, Cp = {N^ - l)/{2Nc) is 
the eigenvalue of the Casimir operator in the fundamental 
representation of the SU(iVc) group, n/ is the number of 
active quarks, (3o = 11 — 2n//3 is the one-loop /3-function 
coefficient, and n > 1. Note that in Ref. |^ the fac- 
tors 3/71 and n/3 were not included in the nondiagonal 
anomalous dimensions Jqq and Jqq presented above, as 
well as a factor 2 was missed in Jqq- The detailed dis- 
cussion of these anomalous dimensions in the one- and 
two- loop approximation can be found in Ref. [pBl . These 



7's allow to define the p parameters entering in Eqs. ( Al) 
and (|A2|) : 



-.71 



(9) ^ 6 r+ iQQ (,) ^ 1 iGQ 

The numerical values of the 7's and p's needed in the 
numerical analysis are given below for QCD with = 3: 



50 
81 ' 

-^2 _ _11 
7gG — 9 ' 

~ -48 
'+ — 81 ' 

P2 — — ' 



-^2 ^ W_ 
IGQ — 243' 

7qg — 3' 

2 101 

'- — 81 ' 



(AlO) 



90' 



The parameters Bn'^ and Bn'' are not determined by 
perturbative QCD and are treated as free parameters, to 
be fixed by data, for example, from the 77' — 7 transition. 



APPENDIX B: THE FUNCTION J(a;,r7) 

The integral J(a;, ry) used in the main text is defined as 
follows: 



J(w,77) = 



dx 



1 -f Lo{x — x) — 2xxr) + ie 



(Bl) 



It is easy to see that this function is symmetric on its 
first argument: J(— r;) = 3(uj,r]). 

The result of the integration depends on the correlation 
between uj and 77 and is the following: 



22 



J(w,?7) 



-(l-r,)2 



arctan 
In 



_2+t£_ 



arctan ■ 



7] — UJ 



^1_^2_(1_^)2 
— iTT 8(77) 



+ (1 -77)^ < 1; 

U^ + il-T])^ > 1, 



(B2) 



where 6(77) is the unit step function. In the hmit uj ~* ±1 this function has a logarithmic divergence 



J(l-2£,77) = 



1 



2(1-7;) 



1 



2e 



(1 - 77)2 



In 



(1-^)' g »^e(77) 

£ 2(1-77)3 2(1-77) 



1 



2£ 



Another asymptotics of the function J (a;, 77) is also used: 



1 



1 - 2£ ' 1 - 2£ 



2(1-7,) 



2£(1 + 773)"1 



In 



(1-7^)2 £(1-277) 77re(77) 



(1-77)2 2(1-77) 



(1-77)2 
_ 2£(l+77) 

(1-77)2 



+ 0{s^). (B3) 



-0(£2). (B4) 



At small values of the second argument 77 — > this function has the following asymptotic expansion: 



J(W,77) ~ 

1 

^ 2^ 



377 7,2(15-4^2) 5773(21-11^2) 



In 



2^2 



l-LO 



— i'!TQ{ri) 



24^6 

77 ^ 772(3 - t^2) ^ 7y3(5 _ 3^2) ^ ^4(35 „ 3Q^2 ^ 3^4) 



(B5) 



2w4 



2tj6 
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